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Dipole-dipole interactions in a square planar array of sub-micron magnetic disks (magnetic dots) 
have been studied theoretically. Under a normal magnetic field the ground-state of the array un- 
dergoes many structural transitions between the limiting chessboard antiferromagnetic state at zero 
field and the ferromagnet at a threshold field. At intermediate fields, numerous ferrimagnetic states 
having mean magnetic moments between zero and that of the ferromagnetic state are favorable en- 
ergetically. The structures and energies of a selection of states are calculated and plotted, as are the 
fields required to optimally reverse the magnetic moment of a single dot within them. Approximate 
formulae for the dipolar energy and anhysteretic magnetization curve are presented. 



I. INTRODUCTION 

Recently some peculiar properties of sub-micron mag- 
netic particles (magnetic dots) fabricated from such soft 
magnetic materials as permalloy, Co, etc., and forming an 
artificial lattice have attracted great attention, see Refs. 

These magnetic dot arrays constitute promising ma- 
terial for high-density magnetic storage media. The dis- 
tribution of magnetization within the dots is quite non- 
trivial. In the absence of an external magnetic field, a 
small enough non-ellipsoidal dot exhibits a single-domain 
nearly uniform magnetization state, either a so called 
flower state or a leaf state.Q On increasing the sizei-pf the 
dot above a critical value, a vortex state occurs.0 This 
vortex state has been experimentally observed ( see Refs. 
^, |5|, H - |l0|) for circular disk-shaped magnetic dots with 
diameters 2R = 200 - 800 nm and thickness L = 20 - 60 
nm. In Rcf. [If], magnetization reversal for an array of 
disk-shaped dots under the influence of a magnetic field 
applied in the plane of the dots has been investigated 
experimentally. In this planar geometry the main con- 
tribution to the total magnetization comes from the in- 
ternal reorganization of each dot's magnetic structure, in 
particular, by displacement of the vortex from the dot 
center, leading eventually to annihilation of the vortex 
at the rim of the dot. In this process, dipolar interaction 
between the dots does not play an essential role. 

In the present work, another case, namely, the ground 
state of an unbounded planar square lattice of thin cir- 
cular disk-shaped dots in an external magnetic field per- 
pendicular to the plane of the dots will be considered. 
We will show that the situation in this case is very differ- 
ent from that in which the field is applied in-plane: the 
main contribution to the total magnetization is now de- 
termined as much by the dipolar interactions of the dots 
as by the external field. In a perpendicular field, the 



dipole-dipole interaction between the dots results in a 
complex specific phase diagram. In particular, a cascade 
of phases with different patterns of dot magnetization 
has been found; these constitute the sequence of ground 
states as a function of the external magnetic field. 



II. MODEL 

In order to formulate the model we need to discuss 
briefly the character of M, the magnetization distribu- 
tion within a single dot in the vortex state. In circu- 
lar cylindrical coordinates, M = M(z,r, x). In suffi- 
ciently thin material, such as that being considered in 
the present work, M does not depend significantly on 
the z-coordinate along the normal to the dot, so that 
we may more simply write M — M(r,x), where r and 
X are the polar coordinates in the dot plane. Then 
the Cartesian components of M for the vortex state in- 
side the dot, M x = M s sin 6* cos (,9, M y = M s s'm9 sm ip, 
M z = M s cos 0, where M s is the saturation magnetiza- 
tion, are determined by ansatz 

6 = 6(r),<p = X + <Po- (1) 

Such a distribution is typical for magnetic vortices 
with a topological charge (vorticity) equal to one, in two- 
dimensional easy-plane ferromagnets. The function 9(r) 
is determined by an ordinary differential equation and its 
solution can be easily found numerically - see the general 
discussion of magnetic vortices in Refs. [l2L O . 

For the theory of magnetic dots made of soft magnetic 
materials, the crystallographic easy-plane anisotropy is 
negligible, and the demagnetizing field H m plays the 
main role. The sources of H m are both the volume 
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"magnetic charges", proportional to divM, and the dis- 
continuity in the normal component of magnetization 
at the surface of the sample (the surface "magnetic 
charges"). The vortex distribution ([l]) has an advan- 
tage compared with others, because for (|l|) divM = 

M s costpo ( — cos H — sin6M, and divM = at ipo = 
\dr r J 

±7r/2. Accordingly, the volume magnetic charges vanish 
and the sole source of the field is the z— projection of H m 
onto the faces of the dot. The states with <p = +ir/2 
and ipo = — 7r/2 have the same energy, i.e. the vortex 
state of the dot is twofold degenerate with respect to the 
sense of the magnetization rotation. In the case of a thin 
enough dot, this gives H m = —4irM s z, i.e. there is an 
effective easy plane anisotropy w m = 27rM 2 cos 2 0. 

The function Of) can be obtained by applying well- 
known methods for treating magnetic vortices in easy 
plane magnets. In the center of the dot (r = 0) the func- 
tion Of) is restricted to two possible values 0(r = Q) = 
0, 7r ; so that cos = p = ±1. Here p is the so-called vor- 
tex polarization (the second topological charge) £j The 
characteristic scale of the variation of the function Of) 
coincides with the value of the exchange length Ao 

A = y/A/A-KMj, (2) 

where A is the inhomogeneous exchange constant. 

For r ^> Ao in the vortex solution, the value Of) tends 
to 7r/2 exponentially. If the dot radius R 3> Ao, then one 
can obtain an acceptable solution in which the limiting 
value Of) = 7r/2 is reached at r = R, the rim of the dot. 
In this case w ir/2 within the major part of the vol- 
ume of the dot, where Aq < r < R. Such states have 
been discussed in a theoretical treatment of magnetic 
dots that included an exact treatment of the magnetic 
dipole-dipole interaction.!;] That treatment showed that 
the out-of-plane magnetization M z is significantly differ- 
ent from zero only in the core region, r < A , and that 
the total magnetic moment of the dot p, has one of two 
values: 

ft = p^z,p = ±l,fj, = 2Tv^LAlM s , (3) 

where £ is a multiplicative constant of the order of 1; in 
fact £ -> 1.361 as A /R -> 0. 

Thus, we arrive at the following simplified picture of 
the state of the single dot. In the greater part of the dot, 
the magnetization lies in the dot plane, rotating about 
the center of the dot, and so, because of its circular sym- 
metry, does not contribute to the total moment \x of the 
dot 13 The state of the dot is fourfold degenerate, with 
ipo = ±7r/2 and the core polarization p = ±1. The value 
of ipo does not manifest itself in the magnetic moment 
of the dot and, as a consequence, has no influence on 
dot interactions. The magnetic moment is directed per- 
pendicular to the plane of the dot and has either of the 
values ±/iz. 



With regard to dot interactions, all dots are in one or 
other of two states: "up" and "down". Since the core 
volume is much smaller than the dot volume, the core 
magnetic moment is small compared with the saturation 
moment of the dot. Although-the dipolar interaction be- 
tween dots is not very strongjLj it is nevertheless the sole 
source of interaction within the dot system. Because both 
the dipolar dot interaction field and the external mag- 
netic field are very much lower than the effective fields 
(exchange and demagnetizing) internal to the dot, one 
can regard the magnetization distribution inside the dot 
as practically unaffected by them. 

An important consequence of this robustness of the 
vortex state is that the polarization p and moment \x of 
a dot remain unchanged under the application of_small 
enough external magnetic fields H e parallel to zS3 It is 
clear that not only the state with H e \\p (light vortex) 
is stable, but also the state H e with antiparallcl to z 
(heavy vortexL remains constant and metastable up to 
\H e \=4wM s m 

Following from the above discussion, the Hamiltonian 
of a dot array can be presented as follows: 

where [i is the moment of a single dot, Pr= ±1, Z, /' are 

dot positions in a square lattice, / = a fix + ny), m,n = 
0, ±1, ±2, ... are integers, a is the interdot distance, and 
H is the external magnetic field parallel to z. The first 
term describes the dipole-dipole interaction of the lattice 
of magnetic dots, the second term is the Zeeman energy. 

It should, perhaps, be stressed that the system mod- 
elled here, consisting of a square lattice of discrete 
dipoles, normal to the lattice plane, with only dipolar 
interactions, is very different from a continuous thin film 
with perpendicular anisotropy. Numerical treatments of 
such films, when performed using a square lattice dis- 
cretization, eg. see Refs. [n| [l?], give rise to a Hamil- 
tonian that bears a superficial resemblance to that of 
the present system, but the essential continuity of the 
film, expressed in the exchange coupling between near- 
est neighbour elements of the numerical discretization - 
the dominant interaction in any sufficiently refined dis- 
cretization - necessarily results in magnetization patterns 
(stripe domain structures) that are entirely different from 
the patterns of discrete moments arising from pure dipo- 
lar interactions between discrete dots on a square lattice 
that are reported here. 



III. MAGNETIC GROUND STATES 

The dipole-dipole interaction is long-range, and it is 
not obvious a priori what structure will constitute the 
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ground state. For a system of particles with dipole- 
dipole interactions and in zero applied field, a theorem 
states that in the ground state the overall magnetic mo- 
ment is zero, and this results in a specific antiferromag- 
netic (AFM) state.ta For instance, in the case of a three- 
dimensional simple cubic lattice of spherical particles, 
a four-sub-lattice structure with non-collinear magnetic 
moments is optimahEj This is possible, however, only if 
the magnetic moment of each particle is free to point 
in any direction. Here, because of the robustness of the 
vortex state supporting the dot magnetic moment, we 
have p = ±1 and uniaxial, giving rise to a quasi-Ising 
model. Determination of the ground state is reduced to 
geometrical considerations. 

It is convenient to divide the initial simple square lat- 
tice, on which the dots are located, into elementary mag- 
netic cells of rectangular shape with (k x I) dots, so that 
the overall spatial arrangement of up and down [p = ±1) 
dots (what we call their "structure", "configuration" or 
"pattern") can be produced from a single such cell by 
a translation T = a(Nx + My), where M — km and 
N = In and m,n = 0, ±1,±2, ... are integers. This 
is appropriate for a magnetic structure with sublattice 
number (i.e. smallest rectangular unit cell size) less than 
or equal to k x I. Note that as the choice of specific values 
for k and I restricts the range of structures that can be 
represented in this way, the search for the ground state 
must, in principle, extend to all integer k, I. The Zeeman 
energy, however, depends only on the relative numbers of 
dots with p = +1 and p = — 1 in a magnetic cell, i.e. on 
the mean moment per dot (//) (= and, of course, 

on the applied field H . 

We have investigated a substantial number of states, 
namely, all states with k x I = 2, 3 and 4, and many 
states with larger values of k and /. Because the dipolar 
interaction falls off with the inverse cube of the moment 
separation, i.e. quite rapidly, it is clear that distributions 
in which dots of like orientation are well apart from each 
other, while those of opposite sign are as close as possible, 
will be energetically most favourable. In particular, one 
notes that the ratio of the interaction energy of a near- 
est neighbour dot pair, to that of a pair of next-nearest 
neighbours, is V8 : 1. It seems extremely improbable 
that two dots of the minority population (which we will 
consistently take to be the down dots with p = — 1) could 
ever be nearest neighbours in a ground-state configura- 
tion. In zero magnetic field, the most energetically favor- 
able distribution is the simple chessboard AFM structure, 
see Fig. [|. In this structure neither of the two (equal) 
populations contains any nearest neighbour pairs of par- 
allel dots. Other AFM structures in which dots with the 
same value of p, do occur as nearest neighbours possess 
very much higher energy as is illustrated by the examples 
in Fig. |. 

For the AFM structure, the energy does not depend on 
the applied magnetic field, whereas for the ferromagnetic 
(FM) structure (with p, = +pz for all dots) this depen- 
dence is maximal. The total mean energy per dot can 



be written as W = W m — (p) H, where W m is the mean 
energy, per dot, of the dipole-dipole interaction and (p) 
is the overall average moment per dot of the distribution 
(zero in the AFM case). In virtue of this, the energies 
Wfm and Wafm of the FM and AFM states, for which 
W m = Wfm, Wafm respectively, become equal at some 
field H* = 2(Wfm — Wafm)/p, and, if these were the 
only states possible, one could expect a first order phase 
transition from the AFM to the FM structure at H = H * . 
However, numerous other states are possible, and the sit- 
uation is very much more complicated. 

In the intermediate region between AFM and FM, 
numerous more complex structures with < (p) < \i 
(" ferrimagnetic" structures) may occur. For these states, 
the dipole-dipole interaction energy, W m , is higher than 
Wafm, that for the optimal chess AFM, but the total 
energy W m .H = W m — (p) H is reduced with increasing 
field H. Therefore, such structures may constitute the 
ground state at finite magnetic fields. To describe these 
structures, it is convenient to introduce the dimensionlcss 
magnetization m = (/i) /p. 

In order to determine the values of H that fix the lower 
and upper bounds of such ferrimagnetic ground states, 
we have calculated the change in dipole-dipole interac- 
tion energy AW m that occurs when the magnetic mo- 
ment of a single dot is reversed in the FM and chessboard 
AFM structures. A simple analysis shows that this en- 
ergy change is determined by the energy per dot in the 
initial states. Reversing the magnetic moment of one dot 
in the FM state requires AW m = AWfm, and in the 
chessboard AFM state AW m = 4\Wafm\- The value of 
Wfm found numerically is 4.516811/x 2 /a 3 . Including the 
magnetic field, one can show that the total energy Wfm 
of the pure FM state and that of the same state, but with 
one dot reversed, coincide at the value H = Hi, where 

Hi = 2W F m/v = 9.033622 ///a 3 . (5) 

Evidently, for H > Hi the FM structure is the most 
favorable, but with H < H±, some magnetic moments 
tend to reverse. When H < Hi , but close to Hi, the 
density of these reversed moments will be very low and 
(m) ~ M- 

The chessboard AFM state can be treated in the same 
way. One obtains Wafm — — 1.322943/i 2 /a 3 and for 
the corresponding threshold field Ho, above which it be- 
comes favorable to switch a dot from antiparallel to par- 
allel to the field direction, one finds H Q = 2\Wafm\/ ^ — 
2.645886 fi/a 3 . 

Hence, intermediate phases with < (/i) < /i exist 
within the finite range of fields, Hq < H < H-y, and 
(fi) — * as H — > H , and (p.) — > 1 as H — > H\. Let 
us consider the nature of the ground states in this field 
range. In the limit of low fields, these states are ob- 
tained from the chessboard AFM by reversing the mag- 
netic moments of a small fraction of the down dot popu- 
lation, leaving the remainder undisturbed. At high fields 
H w Hi, the initial structure is the FM one. In both 
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cases, one expects the flipped dots to be dispersed as far 
from each other as possible, in order to minimize their 
contribution to the dipole-dipole interaction energy of 
the system. If it were not for the constraints imposed 
by the square lattice on which all the dots are located, 
one would therefore expect these flipped dots to lie on an 
equilateral triangular lattice. 

Guided by these considerations, we have sought and 
found excellent candidates for those configurations for 
which W m is minimal, for a number of fixed values of 
m = (fi) /// between and 1. The elementary rectan- 
gular magnetic cells that represent a selection of these 
" optimal" (we drop the parentheses hereafter) configura- 
tions are depicted in Figs. |^ and |^, together with their 
corresponding values of m and W m (with W m expressed 
in units of Li 2 /a 3 ). These W m (also the values for Wfm 
and Wafm given above) were evaluated numerically by 
summing the contributions to the field, at each dot in 
the cell, of all dots within a radius 10, 000.5a. The con- 
tribution from all more distant dots was approximated 
by attributing a uniform areal dipolar moment density 
mix/ a 2 to the area of the plane of the dots outside that 
radius. All values of W m obtained in this way are plotted 
as the points in Fig. ||. 

The sequence of configurations in Figs. [2] and || repre- 
sent some of the (very many) stages in the anhysteretic 
magnetization of the dot array from the demagnetized 
AFM to the fully magnetized FM state. Owing to the 
stability of the vortex state in the dots, it is evident 
that these states probably cannot be accessed sequen- 
tially merely by increasing the applied field. They rep- 
resent stages in an ideal anhysteretic sequence, probably 
accessible only by thermal, or quasi-thermal (e.g. mag- 
netic " shaking" ) cooling through the Curie temperature 
(or quasi-Curie temperature), under the appropriate con- 
stant normal magnetic field H(m). 

It is also evident that the configurations reported here 
constitute only a small sample drawn from an infinite 
sequence of such optimal configurations over the range 
< m < 1: for every rational value of m in this range, 
there exist, in principle, numerous different configura- 
tions, one (or possibly more) of which must possess the 
lowest value of W m . (Henceforward, unless the contrary 
is explicit, W m will be used exclusively to refer to this 
lowest energy for given m.) 

Less evident, but at least extremely plausible, is the hy- 
pothesis that W m increases monotonically with m. Con- 
sider the optimal configurational state, in zero applied 
field, for any specific reduced magnetization m. In both 
the majority "up" dot and minority "down" dot popula- 
tions, the dots are occupying the " energetically best" lo- 
cations available to them. However, by virtue of being in 
the minority population, even the least favorably located 
of the " down" dots is surely more favorably located than 
the least favorably located "up" dot: it has more dots of 
opposite sign with which to interact, and need have no 
dot of the same sign for a nearest neighbour; the least fa- 
vorably located "up" dot, by contrast, is certain to have 



another "up" dot alongside. How, then can it be energet- 
ically favorable to reverse the moment of a "down" dot, 
thereby creating yet another "up" dot? Indeed, this argu- 
ment can be pushed a stage further: not only must W m 
increase monotonically with m, but so must its rate of 
increase, dW m /dm, because the dots being reversed are 
selected in order of increasing stability and new sites for 
reversed dots are increasingly less favorable. It appears, 
however, that dW m /dm, though monotonically increas- 
ing, is discontinuous. For example, consider a state with 
a simple structure like that for m = 1/2 in Fig. |^. It 
is evident that the energy increase on reversing one of 
the minority down dots is substantially greater than the 
energy decrease on reversing one of the majority up dots 
(taking into account the decrease in the former energy 
change and increase in the latter on optimizing the two 
new states). Indeed we have carried out this procedure 
for all but one (that for m = 27/28 which is very close to 
the FM state) of our optimal configuration candidates, all 
of which support this prediction. This aspect is further 
discussed in Section IV. 



IV. ANALYTIC ASYMPTOTIC 
APPROXIMATIONS 

Analytic formulae designed to approximate the dipo- 
lar energy W m in the limits m — > and m — > 00 will 
now be derived. These formulae are both instructive and 
in remarkably good accord with the numerical values of 
W m calculated for specific states and represented by the 
points plotted in Fig. ||. 



A. Approximations near FM state 

1. Equilateral triangular superlattice 

Consider states with m = 1— e, < e < 1. In this 
region, as mentioned earlier, one expects the minority 
down dots to be distributed as far from each other as 
possible, thereby minimizing their positive dipolar inter- 
actions with one another and maximizing their negative 
interactions with the majority up dots. For a given den- 
sity of minority dots, as prescribed by e, their maximum 
separation is known to be ideally accomplished when 
those dots form an equilateral triangular lattice (ETL 
hereafter). In the present case this cannot be precisely 
achieved because all the dots are constrained to lie only 
at points on the fundamental square lattice of spacing a. 
However, when e is very small, the spacing of the minor- 
ity dots A>o, and they can adopt a fair approximation 
to an equilateral triangular distribution, and indeed, as 
e — * 0, this approximation will become very good. In 
developing our approximate formula, we will therefore 
assume the minority dots all lie on an ETL with spacing 
A(e). 
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In order to proceed, we need to know the interaction 
field and energy for dots in an FM state on an ETL. This 
was calculated numerically, in a manner similar to that 
used for distributions on a square lattice described above, 
again summing over a circular region, of radius R = 
10000. 5A, surrounding a central dot at which the field 
of the others is calculated. At lattice spacing A, there 
are (2/-\/3)/A 2 dots per unit area. The region outside R 
was represented, as before, as uniformly polarized with 
the mean dipole moment per unit area, (2 / \f?>) /J, / A 2 . The 
numerical calculation yields an energy per dot Wfma = 
5.517088 [i 2 /a 3 . For A = a, Wfma is substantially higher 
than the value Wfm = 4.516811 /j, 2 /a 3 obtained for the 
square lattice, but that is because the dot areal density 
is a factor 2/V3 higher. For the same dot areal density, 
1/a 2 , we require A = [2/V3] 1/2 a. Because of the inverse 
cube interaction law, the energy per dot, at the same dot 
density, 1/a 2 , is a factor [(V3)/2] 3 / 2 lower than 5.517088 
H 2 /a 3 . This gives Wfma = 4.446373/^ 2 /a 3 for an ETL 
of moment density \ij a 2 , about 1.56% lower than Wfm 
for the square lattice, demonstrating the small but sig- 
nificant energetic advantage of the former configuration. 
(Indeed, the smallness of this energy difference, for such 
very different configurations, is reassuring, for it indi- 
cates that the small departures from the ideal ETL, that 
are imposed by conformity with the underlying square 
lattice, will not introduce any substantial error.) Analo- 
gous to Hi = 2Wfm/ H = 9.033622 /i/a 3 , we will write 
H A = 2W F ma/v = 8.8927451 /i/a 3 . The self-energy per 
dot, of an ETL of dots, with moment /i per dot, and dot 
areal density (e/2a 2 ), is {e/2) 3 / 2 W F MA in the FM state. 
Returning to the m = 1 — e system, we can regard it 
as the superposition, on the uniform square FM dot lat- 
tice, of an (approximately) equilateral triangular system 
of " double- dots" of moment — 2[i and of moment areal 
density — e^/a 2 , i.e. with spacing A = [(2/v / 3)(e/2)] 1 / 2 o 
appropriate to a dot areal density — e/2a 2 . We superpose 
double-dots in order, in effect, to reverse the moments fi 
of that fraction, e/2, of dots on the basic square lattice 
that constitute the (approximate) ETL of up dots that 
require to be reversed to achieve the required overall re- 
duced moment m = 1 — e. The contribution to the energy 
of each double-dot, of moment — 2/i, due to its interaction 
with all the other double-dots, is 4W FM A(e/2) 3/2 . This 
contributes 4(fiH A /2)(e/2) 5 / 2 to the mean dipolar en- 
ergy per dot of the overall square lattice, in this analytic 
asymptotic approximation for m w 1, denoted W a i- In 
addition, the double-dots also experience the field —Hi 
of the underlying square FM lattice on whose dots they 
are superimposed. Consequently the interaction energy 
of the triangular and square lattices is —fieHi, per square 
lattice dot. The self-energy of the square FM lattice is, 
of course, just Wfm = fJ-Hi/2 per dot. Adding the three 
contributions gives for the mean energy per dot of the 
m =1 — 6 system: 



This formula is represented by the curve Wal(m) that 
extends from m = 0.5 to m = 1 in Fig. |J. Agreement 
with the points W that represent the numerical values 
calculated for specific optimal structures is remarkably 
good, not only near m — 1, but over the whole of this 
range. One notes also that over the whole range, the 
approximate values W a \ < W m , the "exact" numerical 
values for specific configurations. This is as it should 
be, because, in the approximation, the minority dots are 
located on an ideal ETL whereas, in the configurations 
treated numerically, the minority dots are restricted to 
points on the square lattice that are close to, but not 
precisely at, the ideal locations. 

2. Square superlattice 

Examination of the specific configurations treated nu- 
merically and illustrated in Fig. |^ reveals that, in some 
cases, notably those for m = 3/5 and m — 4/5, the re- 
quirement that all dots lie on the basic square lattice is so 
restrictive that the minority dots are obliged to lie on a 
square superlattice, no better approximation to the ideal 
ETL being available. It is instructive therefore to mod- 
ify the above treatment and adapt it to a square lattice. 
Whereas the minority dots can only ever approximately 
conform to an ETL, they can lie precisely on a square 
lattice whenever the area per minority dot is k 2 + 1 2 with 
k, I integers. Of the three contributions to the mean en- 
ergy per dot expressed in Eq. ([|), only the self-energy 
of the double-dot lattice requires modification: in place 
of Ha in that equation, we require H\. Or, expressed 
rather in terms of Wfm — 0.5/xifi, we obtain for the en- 
ergy per dot of the system with minority dots on a square 
superlattice: 

Wa = W FM [1 -2e + 4(e/2) 5 / 2 ). (7) 

This expression yields precisely the same values as those 
found directly numerically for the specific structures pro- 
posed for m = 3/5 and m = 4/5 in Fig. || and also, 
for m = 0, that quoted above for Wafm- We have 
Wafm = W FM [{1/{V2) - 1]. 

Because Ha and Hi differ by only 1.56% and the 
double-dot lattice self-energy term is proportional to 
e 5 / 2 , Wo differs very little from W a \ over the range 
1/2 < m < 1; the difference in the worst case, m = 1/2, is 
only 1.11%. We do not include a curve representing Wn 
in Fig. [I| because it can scarcely be distinguished from 
that for W a i- Over the range 1/2 < m < 1, whereas W a i 
represents a close underestimate of W m , W u provides a 
similarly close overestimate. 

B. Approximation near AFM state 



W a i = M #i(0.5 - e) + 4( Mj ff A /2)(e/2) 5 / 2 . (6) 



Here we consider states with positive m close to 0. 
The treatment resembles that in the environs of the FM 
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state discussed above. We again expect the dots that 
depart from the chessboard AFM structure ("up" dots 
this time), to be distributed as far from each other as 
possible, in locations approximating an ETL. Again we 
assume these "exceptional" dots all lie on an ideal ETL 
with spacing A, and replace e in the above discussion 
of the triangular lattice energy by to. We now superpose 
the triangular lattice of " double-dots" on the AFM lattice 
instead of the FM one, an important difference being that 
we must, of course, place the positive double-dots only 
on top of negative dots of the underlying AFM lattice, 
whereas all dots in the FM lattice were equivalent and 
available for reversal. 

Apart from the replacement of e by to, the expression 
for the self-energy of the double-dot triangular lattice 
is the same: 4(fj,H/\/2)(m/2) 5 / 2 per dot of the overall 
square lattice. However, the interaction energy of the 
triangular and square lattices is now positive, +TO/ii?o 
(as against —e/j,Hi) per dot of the square lattice and the 
self-energy of the square AFM lattice is —fj.Ho/2 per dot. 
Adding the three contributions gives for the mean energy 
per dot of the overall lattice of weak reduced magnetiza- 
tion m: 



W aQ = »H (m - 0.5) + 4 m (H a /2)(to/2) 5 / 2 . (8) 

Note that the two asymptotic approximations W a o and 
W a i happen to coincide at to = 0.5. The expression for 
W a o is represented by the curve Wa0(m) that extends 
from to = to to = 0.5 in Fig. Again, agreement 
with the points W that represent the numerical values 
calculated for specific optimal structures is remarkably 
good, not only near to = 0, but over the whole range 
< to < 0.5 . However, the agreement is not quite as 
good as was the case for W a \ , very probably owing to the 
additional constraint that only dots from the down-dot 
population are available for reversal, whereas in the FM 
case, all dots were available. This makes it somewhat 
harder (in an actual configuration, but not, of course, 
in the analytic approximation) to arrange the non-AFM 
dots close to the ideal triangular lattice. This cannot, 
however be the only reason for the lower agreement be- 
cause, whereas for the approximation near the FM state 
all W a i < W m , some W a o > W m . An outstanding exam- 
ple is that for to = 1/3 where W a o — W m > 0.01/Lt 2 /a 3 . 
The reason for this is that, in an actual configuration, 
the constraint, assumed in the analytic treatment for 
to < 0.5, that the ideal configuration is obtainable by re- 
versing only some negative dots in the chessboard AFM 
configuration and without any further rearrangement of 
the structure, does not apply and, for values of m suf- 
ficiently higher than zero, a lower energy than W a o can 
sometimes be achieved by violating this supposed con- 
straint - see, for example, the simple configuration for 
to = 1/3 in Fig. |. 



V. ANHYSTERETIC MAGNETIZATION 
CURVE 

Consider two optimal configurations, one of reduced 
magnetization to, the other of higher magnetization 
(to + dm). Their energies in a normal magnetic field 
H, W m>H = W m - m^iH and W {m+Sm ^ H = W {m+Sm} - 
(jn + Sm)fiH , will be equal only in the field -ff m ,(m+5m) = 
(W(m+8m) ~ Wm)/(p5m). In the limit 5m — > 0, we obtain 
the following differential expression for the anhysteretic 
magnetization curve H m (m): 

fj,H m = dW m /dm. (9) 

(Our remarks above concerning local discontinuities in 
dW m /dm and its monotonic increase with to are clearly 
relevant here also to H m .) 

In Fig. [5] we have plotted points representing approxi- 
mate values of H m (to) derived from the numerical values 
of W m calculated for our set of candidates for optimal 
configurations. They are labelled dW/dm, and repre- 
sented by circles. These are approximations to H m given 
by (W m ( i+1 )-W m (i))/n{m( i+ i) -to (i) ) and plotted at the 
values to = (to(j +1 ) + m^)/2 located midway between 
successive configurations of the known set. (All fields are 
shown in units of [i/a 3 .) 

Also plotted in Fig. ||, where they are labelled Hal and 
HaO, arc smooth curves representing H a \ and H a ( h the 
analytic approximations to H m obtained by differentiat- 
ing the asymptotic equations (^|) and (@) and dividing by 

H al = -^dWax/de = Hx- (5/2)ff A (e/2) 3 / 2 , (10) 

H a0 = n^dWao/dm = H + (5/2)H a (to/2) 3 / 2 . (11) 

(The corresponding equation for H a = —fi dWn/de is 
not plotted in Fig. |5|: it is practically indistinguishable 
from H a i over the appropriate range 0.5 < m < 1.) Both 
Equations (|l^) and ( pd| ) provide a good fit to dW m /dm 
over their appropriate data ranges, < to < 0.5 for H a o 
and 0.5 < to < 1 for H a \. As one would expect for such 
asymptotic approximations, the fits are best towards the 
limits to = and to = 1. Unlike the expressions (|^) 
and (^) , from which they are derived, equations (|lC]) and 
( |i"T| ) do not lead to coincident values at m = 0.5. That 
expressions (||) and (||) should yield the same value for 
W m at to = 0.5 is a remarkable coincidence; that their 
slopes dW m /dm should also agree there, is scarcely to be 
expected! In fact the discontinuity between them is in 
remarkably good agreement with the corresponding step 
in the data for dW m /dm in the vicinity of m = 0.5. 

VI. CONFIGURATIONAL STABILITY 

The fields H + (to) and H- (to) required to render en- 
ergetically favourable the reversal of a down or up dot 
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in the optimal configuration for any value of m, and the 
difference between them, H + — H_ , provide measures 
of the stability of that configuration. Limiting exam- 
ples of these fields, for the AFM and FM states, namely 
Hq = H + (m = 0) and Hi = i?_(m = 1), were discussed 
earlier in this article. (In fact, the AFM state also has 
H-(m = Q) = — Hq, so that it is stable over a wide field 
range, 2Hq.) As mentioned earlier, we have determined 
numerically the extra energies AW rn+ and AW m - re- 
quired to reverse a down dot or an up dot respectively 
(including local rearrangement of the resulting configu- 
rations to minimize their energy) in almost all our can- 
didates for the optimal states of reduced magnetization 
rn. 

The values of the corresponding fields, H±(m) = 
AW m ±/2/j, are plotted in Fig. ||. Some examples illus- 
trating the local re-organization of the dot distributions 
after single-dot reversal are presented in Fig. ^. Interme- 
diate values of m were selected for display in this figure 
because the patterns in that region are rather more com- 
plex and varied than those near to = and m = 1 and 
usually give rise to greater ranges of field stability. 

In Fig. [| it is evident that a salient feature of the fields 
H+ and in this intermediate range of magnetizations 
is indeed the very substantial width of the gap between 
them, H + — H- , that defines the field range of stability 
of the pattern against single dot reversal. Towards the 
limits to = and m = 1, that gap shrinks towards zero. 
However, for the structure of the pure AFM state, pre- 
cisely at m = 0, the gap is at its widest: 2Hq. (_ff_ for 
the pure AFM state is off-scale in Fig. |5| ) Throughout 
the range < to < 1, one expects H + > H m > H-, 
but, of course, we do not have values of H m at the values 
of m that correspond precisely to the optimal configu- 
rations. However, one observes in Fig. || that, for al- 
most every configuration (label i), H + (m^) exceeds the 
next value of H m w {W m{l+1) - W m({) ) / ^i{m [l+1) - m^), 
with m( i+1 ) > 777(j), while, in the reverse direction, 
H—(mm) < H m (i_x\, This means that the field required 
to reverse a single down dot in one of our optimal con- 
figurations, even after allowing for local rearrangement 
of the pattern to minimize the increase in dipole-dipole 
interaction energy, usually exceeds the field required to 
render energetically favourable the reversal of the infinite 
array of down dots required to change (anhysteretically) 
the pattern to that of the optimal configuration appro- 
priate to the next higher value of to. Similarly, to reverse 
a single up dot, also after local rearrangement, usually 
requires a reduction of the field to a value below that 
needed to favour the (anhysteretic) reversal of the in- 
finite array of up dots needed to establish the optimal 
configuration for the next lower value of to. 

Examples may be helpful here. Compare the central 
rectangles for the states to — l/2(+) and to — 5/9 in 
Fig. ^. The former has the lower mean reduced mo- 
ment, 13/24 < 5/9, and is embedded in an infinite ar- 
ray of even lower mean reduced moment, to = 1/2, yet 
it requires a higher field, H + (m = 1/2) = 6.681 > 



H m (m = 0.5277...) = 6.589 (in ^/a 3 units), to equili- 
brate it with the to = 1/2 configuration than does the 
periodic to = 5/9 pattern. Likewise, compare the cen- 
tral rectangles for the states m = 1/2 and to = 5/9(— ) 
in Fig. ^. Both have the same mean reduced moment 
of 1/2, but the to = 5/9(— ) rectangle is embedded in 
an infinite array of higher reduced moment, to = 5/9, 
yet it requires a lower field, H_(m = 5/9) = 5.986 < 
H m (m = 0.5277...) = 6.589 (in fi/a 3 units), to equili- 
brate it with that to = 5/9 configuration than does the 
periodic m = 1/2 pattern. 

The reason for this prima facie somewhat paradoxical 
behaviour is that, when only a single dot is to be switched 
in an optimal regular periodic configuration, only very 
local rearrangement of the dot pattern can reduce the 
energy of the resulting perturbed system whose overall 
pattern, remote from the switched dot, must remain op- 
timised for the original level of mean magnetization. It 
follows that the sequence of points constituting the ideal 
anhysteretic magnetization "curve" must correspond to 
a discontinuous sequence of stable energetically optimal 
dot configurations. Magnetization by monotonic varia- 
tion of the applied field alone will exhibit very substan- 
tial hysteresis and will realize very few, if any, of the ideal 
magnetization patterns. 



VII. SUMMARY AND CONCLUSIONS 

Unbounded two dimensional arrays of thin circular 
disk-shaped magnetic dots on a square lattice have been 
considered in the presence of a magnetic field perpendic- 
ular to the dot plane. The radii and thickness of the dots 
is such that a radially symmetric vortex magnetization 
structure is stable; at radii outside a relatively limited 
core, the magnetization lies almost in plane and parallel 
to the rim. The net moment fi is due to the vortex core 
and is normal to the dot plane; it is practically unaffected 
by normal fields of magnitude comparable to those from 
other dots. 

For every rational value of the reduced magnetization 
of the array, to = (/j,) //i < 1, (ignoring sign) there ex- 
ist very numerous possible arrangements of up and down 
dots, one (or more) of which must have minimum dipo- 
lar interaction energy W m . We have calculated that en- 
ergy for a range of excellent candidates for these ground 
states at various values of to, in particular those of the 
chessboard AFM state of to = and the uniform FM 
state with m = 1. We suggest and argue that, for the 
sequence of optimal configurations, W m and dW m /dm 
increase monotonically with to, and this is supported by 
our data in Fig. |^. (dW m /dm is likely to be locally dis- 
continuous). An analytic formula for W m (m) is derived 
on the assumption that for to — ► 1, the minority down 
dots arc located near points on an equilateral triangular 
super-lattice. A similar formula is derived for to — > 
relating to the excess up dots. Remarkably, these for- 
mulae agree at to = 1/2 (though their gradients differ) 
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and fit the data for specific states very well, especially 
near to = and 1. A similar formula involving a square 
super-lattice instead of the triangular super-lattice gives 
results in precise agreement with the data for specific 
structures with m = 0, 0.6 and 0.8, these being those 
structures in which the minority dots do indeed lie on 
square super-lattices. It differs from the first asymptotic 
formula by little more than 1.11% over the appropriate 
range 0.5 < m < 1. 

Approximate data for the anhysteretic magnetization 
curve H m (m) = ^ 1 (dW m /dm) are derived from the 
data for specific states and compared with the predic- 
tions of the analytic formulae in Fig. The agreement is 
good, particularly near to — and 1. The fields predicted 
by the two asymptotic formulae differ quite sharply at 
to = 1/2, where they indicate a step discontinuity in field 
that matches a similar step in H m (m), the data from the 
series of specific states. 

The stability of the optimal configurations found was 
explored by determining the minimum field H + required 
to reverse a single down dot and likewise the maximum 
field at which a single up dot would reverse, assuming 
in both cases local reorganization of the resulting dot pat- 
tern to minimize its energy. The gap H + — H- between 
the two fields that indicates the stability of the config- 
uration increases irregularly from near zero, for to just 
greater than zero, to a rough maximum around to = 1/2 
and then decreases again towards zero at m = 1. An 
exception to this general trend is the very large value 
at m = where H + = —H- = Ho. For the optimal 
configurations at most values of to, H+ > H m+ , and 
H- < H m —, where H m+ and H m -, refer to the optimal 
configurations studied at neighbouring values of magne- 
tization, just above and just below m. It follows that 
there exists, with increasing m, a sequence of stable opti- 
mal configurations with energy barriers between them. It 
must be stressed, therefore, that the field curve H m (m) 
is strictly an ideal anhysteretic magnetization curve: ow- 
ing to the stability of the individual dot moments, the 



individual dot distributions are likewise very stable, and 
the infinite sequence of energetic optimal states cannot 
be traced experimentally by monotonically increasing (or 
decreasing) the normal field alone. To realize any specific 
state it would be necessary - but perhaps not sufficient - 
to cool the sample through the Curie temperature subject 
to a normal magnetic field of the appropriate strength. 
Another possibility would be to destroy the stability of 
the other metastable phases and achieve the phase with 
minimal energy by a form of magnetic shaking, e.g. by 
the application of fluctuating fields of decreasing ampli- 
tude. 

For magnetic storage applications, on the other hand, 
the stability of the dot configurations is advantageous, 
indeed, essential. Experimental observation of these dot 
arrays and transitions between their states may perhaps 
also be useful for the determination of the basic dot pa- 
rameters, in particular the radius of the vortex core and 
its magnetic moment. 

Last, but not least: these results are clearly not re- 
stricted to a dot lattice of the type considered here, 
but also apply directly to any square lattice of identical 
dipoles that are restricted to the two senses of normal ori- 
entation. They may also be relevant to the description 
of other uniaxial dipole-coupled systems of small parti- 
cles, for example, thin films of granular magnets with 
easy axial anisotropy (shape or crystallographic) with 
a perpendicular easy axis and negligible exchange cou- 
pling between granules. Such properties are characteris- 
tic for thin films prepared by simultaneous evaporation 
of permalloy .ayd silver with small enough concentrations 
of permalloyEJ 
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FIG. 1: Different antiferromagnetic states for a square mag- 
netic dot lattice. The dipolar energies per dot of the lattices 
(in units of (/i 2 /a 3 - see text) have been found numerically 
and printed below each sketch. Open circles o denote the 
"up" dots with moments p = fj>z; solid circles • correspond to 
"down" dots with (1 = —fj,z. 
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FIG. 2: Some lattice structures with intermediate values of 
rn — (fi z ) /[i, < m < 5/9, and thought to have least dipolar 
energy for that value of m. The values of m and dipolar energy 
per dot of each structure (in units of /J? /a 3 ) are printed below 
it. 
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FIG. 3: As Fig.2, but for 3/5 < m < 1. 




FIG. 4: Dependence of optimal lattice dipolar interaction 
energy per dot, W m , on reduced lattice magnetization m. The 
points W represent values calculated for specific structures 
(as in Figs. 2, 3). Curves Wal and WaO represent the analytic 
asymptotic formulae in Eqs.(H), (H). (All energies are in units 
of M 2 /a 3 .) 
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FIG. 5: Anhysteretic magnetization curve for magnetic dot 
lattice. Solid circles labelled dW/dM represent approximate 
values of the equilibrium field H m (m) = n~ 1 dW rn /dm de- 
rived numerically from the dipolar energies of the succes- 
sive optimal lattice structures studied (see text). Threshold 
fields for reversing a single dot, H- and H+, are also shown. 
Curves Hal and HaO follow the analytic asymptotic formulae 
for H m (m) in Eqs.(|Io[ 0). (All fields are in units of n/a 3 .) 
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FIG. 6: Examples of optimal re-ordered structures after a 
single dot has been reversed in an optimal periodic lattice con- 
figuration of specific magnetization m. The central column 
shows the structures of individual cells in the initial periodic 
rectangular superlattice. Reversing a single up dot in a single 
such cell ("central" cell) and then rearranging to minimize 
the energy, alters the central cell structure to that in the left 
column. Similarly, the result of reversing a single down dot 
is shown in the right column. None of the other cells in the 
superlattice is altered. 



